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Abstract. The electromagnetic wave field propagating in a helical wave guide 
is decomposed in an angular momentum basis. Eigenmodes are calculated using 
a truncation in I and a discretisation of the boundary condition. Modes slightly 
slower than light are considered for relativistic particle acceleration. 



1 Introduction 

An electromagnetic wave can accelerate charged particles over long enough dis- 
tance if 1) it has a longitudinal electric component 2) its phase velocity is equal 
or close to the particle velocity, therefore less than c. A free plane wave misses 
the first condition and a wave guided by a metallic straight tube miss the second 
condition. The usual solution to satisfy both conditions is a chain of coupled 
resonating cavities. Another solution is the helically bent tube. Here the "pho- 
ton path" is longer than the corresponding path in a straight tube, so that some 
modes can be made slower than light, while having a non-zero longitudinal elec- 
tric field on the helicity axis. 

A helical wave guide is defined by its axis Oz, its period A = 2ir/q and its 
base curve r(s) = (X(s),Y(s)) which is the intersection of the wall with the 
z = plane. The wall surface is generated by helices H s of parametric equation 

(x,y,z) = (K(qz)r(s),z) (1) 

where the 2x2 matrix IZ(qz) describes the rotation of angle qz about Oz. 

The traveling wave tube, used for generation of radar waves, is a special kind 
of helical wave guide. It is made of a cylindrical metallic tube inside which a 
thick conducting helical wire is fixed. Its base is not simply connected : it is a 
disk with a hole corresponding to the section of the wire at z = 0. Therefore 
the base curve has two disconnected pieces : the circular base of the cylinder 
and, inside it, a small closed curve. The same device, with a superconducting 
helix, has been studied for particle acceleration [TJ [5] . 

A tube where the wall itself possesses a helical structure can be made by 
drawing helical grooves on the internal face of an initially cylindral tube. In 
the case where the groove section is made of rectangular steps, the propagation 
modes can be calculated using an approximate field-matching technique [3], or 
a 3-dimensional computer code [4]. 
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Figure 1: Base curve (b.c.) and sections of the helical wave guide at successive 
values of z. The rotating frame {OX, OY} of the last section is shown. In this 
example the base curve is an ex-centered circle of radius a. The balls (s.p.) 
represent the sutture points. 



In this paper we will consider helical tubes which have smoothly curved 
walls, associated to smooth base curves. To calculate the modes we will use a 
truncated partial wave expansion and enforce the boundary conditions only at 
discrete points of the base curve, called sutture points, the number of sutture 
points being equal to the number of partial waves. For definiteness, we will take 
the base cuve to be a ex-centred circle, as shown in Fig[T] The guide parametres 
are then 

• a — radius of the base circle, 

• e € [0,1] — relative amplitude of the helical deformation. The centre of 
the base circle is at distance ea from Oz. 

• q = 2tt / (helix period). 

The average angle of an helix with Oz {i.e., ? minus the pitch angle) is a = 
arctan(qa). Another significant quantity is the maximum helix angle of the wall 
with Oz, (3 = arctan(e<7a). The case (3 = (q = or e = 0) corresponds to a 
straight cylindrical guide. 

1.1 The scalar field case 

We first study the simpler case of a scalar field $(t, x, y, z), which obeys the 
Klein-Gordon equation (in units c = 1) 

{-d 2 t +d 2 x + d 2 y + d 2 z ) x,y,z) = 0, (2) 

with the boundary condition 

$(i, x, y, z) = on the wall . (3) 
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In a fixed z transverse plane we introduce a "rotating frame " whose basis vectors 
rotate about Oz by the angle qz when z increases. The fixed- and rotating-frame 
coordinates, r = (x, y) and R = (X, Y), are related by x + iy — e lqz (X + iY). 
The corresponding polar coordinates are (r, if) and (r, ip — qz) respectively. We 
look for a solution of the form 

$(t,x,y,z) =e iPz - iut V(X,Y), (4) 

P is the longitudinal pseudo-momentum (in the Bloch wave sense). It is the 
eigenvalue of the operator 

"Phciic = Pz + qJz = -id/d z - iqd/d v (5) 

which is conserved due to the invariance under the helical displacement 

{translation by Az} x {rotation by qAz} . (6) 

For a given lu the Klein-Gordon equation with the boundary condition has a 
finite number of solutions or modes. P and lu are linked by the dispersion relation 
P = P m {uj) or lu = io m {P). The phase and group velocities are w p h = w/P, 
Vgf = doj/dP. 

Since $ is an eigenvector of ([5]) we can decompose it in free cylindrical partial 
waves of angular momenum I and longitudinal momentum pi = P — Iq 

S(t, x, y, z) = e— 4 a l ^ e l{P ~ lq)z Am) (r) , (7) 

/— — OO 

where 

k(l) or h = {u? plf 2 = [u 2 -{P- Iqf] 1/2 (8) 

is the transverse momentum of the I component and ipi t k(r) is a solution of the 
radial equation 

V>" + ip'/r + (k 2 - l 2 /r 2 )ilj = . (9) 

ki is real for 

{P-u)/q<l< (P + w)/g, (10) 

1 /2 

otherwise it is pure imaginary : fe; = iki, with ki = UP — lq) 2 — lu 2 ] 

The solution of ([5]) is proportional to the Bessel function J; (kr) or the mod- 
ified Bessel function Ii{nr) = Ji{inr)/i l . We will fix the normalisation of i/ji k{r) 

by 

lfo,k(r) = J W {kr)/k^ = I w {nr)/K^ . (11) 

With k = in the last two expressions are equivalent. The fcl'l denominator 
allows a finite limit for , 0/,fe( r ) when k — > 0. 

Let us write the boundary condition for t — z — 0. II will be satisfied also at 
all t due to the invariance under time translation and all z due to the invariance 
under helical displacement. For every point r(s) = (X(s),Y(s)) of the base 
curve one has the constraint 

]>>e^ s ) J w [kr(s)]/W =0, (12) 
where (r(s), f(s)) are the polar coordinates of r(s). 
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1.2 The truncation-and-sutture method 

One keeps only a finite set C of partial waves and selects an equal number 
of "sutture points" ri,...r N of the base curve (see FigO}. For instance C = 
['infj'sup] with l sup = i; n f + N — 1. The ai are approximatively determined by 
the system of linear equations 



{ r m l Pn) being the polar coodinates of r„. A nonzero solution exists if the 
determinant est zero: 



This is the dispersion relation in the truncation-and-sutture approximation. The 
following conditions may be imposed on C and the r„'s : 

• the r„'s should give a precise enough definition of the base curve, 

• I (ifn+i — fn) should not be too large compared to unity, 

• ki (r n+ i — r„) should not be too large compared to unity. 

The ultimate criterium of the quality of the approximation is the stability of 
the roots of (fT4"|) when N is changed by one or two units. The coefficients a; 
should be stable also, at least for I far enough from the border of C. 

Case of a symmetric base. If the base is symmetrical about Oz 1 the modes 
are either even or odd under this symmetry. One can study the even and odd 
modes separately by taking either even or odd Vs in C. Similarly, if the base has 
a ternary symmetry (invariance under rotation by 2ir/3), C can be restricted to 
either / = or I — 1 or I = 2 modulo 3. To avoid redundant equations in (fT4"|) , 
the points f n should be all different modulo the symmetry. One can choose 
them on one part of the base curve only, the other parts being deduced by the 
symmetry. 

1.3 Numerical example in the scalar case 

We took a circular base with parameters e = 0.4, q = 1 and imposed the phase 
velocity v p h — lo/P = 1 — 10~ 5 . Plotting In |det{C„;}| as a function of wa for 
different sets C and equidistant r„'s, with N ranging from 5 to 21, we found zeros 
at uja — 5.0, 7.0 which are stable at the 0.1 precision when C contains at least the 
set [—1, +5]. Another root, uua — 8.5, seems to stabilise for C D [—1, +7]. Higher 
modes are seen, but would require more partial waves to get stable. Fig. 2 shows 
the scalar field *(A, Y) = $(0, X, Y, 0) of the tua = 5.0 mode. Due to the low 
number of sutture points, the boundary condition ^(X, Y) = at the boundary 
is only appoximately satisfied. One observes that ^(X, —Y) = $>*(X, Y). This 




(13) 



det{C ni } = 0, where C nl = e a *'" J w (kr n )/k^ . 



(14) 
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Figure 2: 3-dimcnsional representation of the scalar field ^(x, y) — <&(0, x, y, 0) 
of the uia = 5.0 mode, for a circular base with e = 0.4 and q = 1. Left: real 
part; right: imaginary part. 



property comes from the symmetry of the base about the X axis. The maximum 
of 1 3?^ | is not at the centre of the base, but farther from the helix axis. This 
is a manifestation of the centrifugal effect. One can also observe that the phase 
arg<3> = arctan(S'I//sftvl/) increases with the azimuth, which means that the 
average angular momentum is positive, as expected. 



2 The Maxwell field 

We replace the Klein-Gordon equation by the Maxwell equations in vacuum, 
V x E = -d t B V • B = 

VxB = 9 t E V-E = 0. (15) 

These equations are invariant under the duality B ^ E, E ^ — B. In the 
following we gather the fields E and B in a 6-component vector F = {E, B}. In 
analogy with (0J we look for solutions of the form 

F(i, x, y, z) = e lPz ~ 1 " 1 K(qz)F(X, Y) , (16) 

in the complex notation (the physical field is the real part of it). T = {£, B} is 
the field in the z — plane and lZ(qz) is the rotation matrix of (TT]) applied to 
the vectors £ and B : 



E x = £x cos qz — £y sinqz , E y = £x sniqz + £y cos qz (idem for B). (17) 

To write the boundary conditions we introduce the vectors 

'X'(s)\ f-qY( S )' 
a(s) = r'(s) = I Y'(s) , f(s) = z + qz x r(s) = qX(s) I, (18) 
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which are tangent to the wall at z = 0, and the normal vector 
v(s) — <j(s) x r(s) = 



Y'(s) 
-X'(s) 
q{XX' + YY') 



(19) 



The boundary conditions are £ ■ a = 0, £ ■ f = and B ■ v = 0. In cylindrical 
coordinates, they read 



£r (J r ~\~ £tp O '^p — . 

qr£ v + £ z = . 
B r a v + (qrB z - B v )a r = . 



(20) 
(21) 
(22) 



These 3 conditions are not independent. Only two suffice. They break the 
duality B -> E, E -> -B. 

F can be decomposed in a manner similar to ([7]). However, since the photon 
has two polarisation states, there are two free partial waves for each I, which is 
now the total (spin + orbital) angular momentum. 

2.1 Study of a free partial wave 

Let us first study separately one partial wave of definite u>, I and p (here p = 
pi = P — Iq). In cylindrical coordinates it takes the form 



( K \ 

E v 

E z 
B r 



e ~iut e ipz e ilip 



( £r(r) \ 

B r (r) 
. B v (r) 
\B z (r)J 



(23) 



In this subsection the quantities written in caligraphic depend only on r. The 
Maxwell equations become 



(l/r)£ z -p£ v = 
—d r £ z + ip£ r = 
d r ) £ v - i(l/r) £ r = 



LU B r 

iui B, 
iui B, 



(24) 



and the dual equations obtained by B — > £, £ 



-B. These equations are 



redundant and one can omit the third one and its dual. 



The TM/TE basis. A first basis is made of the so-called transverse magnetic 
(TM) and transverse electric (TE) states. In cylindrical coordinates, 

/ +ipkJ' w (kr) \ 
-pl/r J\i\(kr) 
= k -\l\- 2Slt0 +k 2 J w (kr) ^ / 25 v 



(r) = 



f£ r (r) \ 

B r (r) 
\B z {r)J 



(TM) 



+k 2 J m (kr) 
+u>l/r J|/| (kr) 
+iu> k J.'j, (kr) 
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with k = (uj 2 - p 2 ) 1 / 2 . The TM and TE fields are dual, up to a factor i 
B( TE ) = -i£(™\ £< TE ) = For the I = waves, 



■0 being defined in (fTTj) . 





ipo 


— icj "01 

V o J 



( °, \ 

+0J1P1 



-pip! 



V -iipo I 



(26) 



The chiral basis. When k — > with Z 7^ 0, jFC™) and ^"( TE ) become equal 
up to a sign and they no more form a complete basis. We avoid this degeneracy 
with the combinations T+ = ^ TE ) + T- = T^) _ j(tm) j which 

are respectively the states of right and left chirality, or helicity \ = +1 and 
X = — 1. They are such that E = ?xB. For real fe a right-handed (resp. left- 
handed) solution is a superposition of plane waves of positive (resp. negative) 
helicity. More precisely, we normalise these states as : 



T 

J L 



XV + siuj 



( -i [ipx-i + (p - x'^) 2 ^a+i] \ 

si [ipx-i - - x'u) 2 ipx+i] 

2(p - x'^) -0a 

-X [V'A-l + (P ~ X'u) 2 1p\+!\ 

-ix' [i>\-i - (p - x'u) 2 ip\+x] 
\ -2«X (p ~ X'w) 0a / 



(27) 



Here si = sign(l), \' = s/x, A = |Z| and 0a 
recurrence relations 



J\{kr)/k x . Use was made of the 



(I fx) Ji(x) = LVxOr) + J,+i(ar)]/2 , J/(.x) = [J^ x (x) - J l+1 {x)}/2 . 

The denominator (xP + Si"-"- 1 ) eliminates a kincmatical zero, so that these states 
have finite and linearly independent limits when p — > — x'w. We will take the 
chiral basis except for the Z = states, where we keep the TM-TE basis (|26p . 



2.2 The partial wave expansion 

According to (|26H27[) . the decomposition ([7]) is replaced by 



F(t,r,<p,z) = e~ luJt e lPz 



(TE) 



+ e-*-' ^ e<*{0« [a, f u , p(1) ,,,^ =+1 (r) + b t ^ mx >=-i(r)] , (28) 

with p(l) = P — Iq. For Z 7^ the a; (resp. 6;) coefficients multiply the states 
where the chirality has the same sign as Z (resp. opposite sign to I). 
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The boundary conditions (|20ti21[) read 

= (p<j r a + iLoa v b ) i/jx 
+ } e tltp { ai[(a r +isia v> )ipx-i + {p - uj) 2 (o- r - isia^ipx+i] 

+ b t [ (oy + isia v ) tpx-i + {p + u) 2 (oy - isicr^,) ip\ + i] } , (29) 

= a V>o + uqr ipt b 

+ bi [siqr (ipx-i - (p + u) 2 + 2(p + w)Va] }• (30) 

1 /2 

In these equations, A = |Z|, p = pi = P — Iq, k = k\ = [uj 2 — pf\ , the 
argument of ip\-i, and ipx+i is h r ( s )-> f — an( A & = &( s ) is given by 

m. 

Generalising the truncation-and-sutture method of Section 1.2 we write con- 
ditions (|29H30[) for N points r„ of polar coordinates (r n ,ip n ), keeping only N 
values of I. It gives 2N linear equations linking N coefficients ai and N coeffi- 
cients 6/. The dispersion relation is the vanishing of the 2N x 2N determinant 
of this system, the coefficients of which are given in Appendix (|35H36p . £ and 
ri, ...rjy are chosen like in the scalar case. 

We have tested this method with a straight (q = 0) cylindrical wave guide. 
The well-known TM and TE lower modes were found using no more than ten 
points, even when the Oz axis was not the guide axis, but shifted by about half 
a radius. 



3 Application to the acceleration of relativistic 
particles 

A particle moving along the axis of the helical wave guide is accelerated or de- 



celerated by the E z component at r = 0, whose physical value is 5R{ao 



APz — iujt 



}• 



Besides, the I = ±1 partial waves give a transverse field at r 
nents of which are, in the rotating frame 

9f{(ai + 6i + a_i + 6_i)e iPz - 
5i{(ai +h- a_! - 6_i) e lPz - lult } 
»{(-ai + h + o_i - e lPz ^ *} 
3{(ai - h + a-i - fe_i) e lPz ^ 1 } . 



0, the compo- 



E x 
E Y 
B x 
B Y 



■iujt \ 
-iujt T 



(31) 



If an electron moves along the axis at the same velocity as the wave (z — v(t—to) 
with v = u p h)i it will experience a Lorentz force f = (— e) (E + v x B) which is 
constant in the rotating frame. In the ultrarclativistic case (u p h ~ 1), this force 
is 



fx 
fr 



(-e) x 25R 



-iujto 



—i 
l 



b-i 



(32) 
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In the laboratory frame, it is the force excerted by a helical undulator having 
the same period as the helix. We therefore get a parasitic undulator radiation. 
It can be avoided with a base which has a symmetry of order S > 2 about Oz, 
choosing a mode involving only the waves / = modulo S. 

4 Numerical results 

Like in the scalar case (Section II -3[) we took a circular base with parametres 
£ = 0.4, q = 1 and imposed the phase velocity v p h = 1 — 10~ 3 . Plotting 
In |det{C„z} or, better, arg |det{C n ;}| we located modes at ua — 2.28, 3.6, 4.8 
and 5.0 which are stable at least at the 0.1 precision for N ~ 15 (they were 
the same for C = [—6, +8], [—5, +9] and [—3, +11]). The first three modes 
seem to be mostly of the TE type : \E z (0)/B z (0)\ = 0.3, 0.25 and 0.4. We 
give below more details about the fourth one (ua = 5.0) which is mainly TM : 
|£,(0)/B,(0)|=5.4. 

Fig. 3 shows the cross section of the fields along the segment x E [—0.6, +1.4] , 
y = 0, which is a diametre of the base circle (in units a = 1). At the extremities 
the boundary condition B r = 0, which was not directly imposed by (|29II30|) is 
well satisfied. 

Fig [4] shows the 3-dimensional plot of \S r (x,y)\, \£ v {x,y)\, \£ z {x,y)\ and 
\B z (x,y)\. We did not plot \B r \ and \B V \, which look like \£ v \ and \£ r \. 

A figure of merit of the helical wave tube as an accelerating device is the 
ratio of the accelerating field over the strongest electric field at the boundary, 

f.o.m. = £ 2 (0)/max(|£(r(s))|) . (33) 

Indeed a too strong field at the boundary may trigger an electron emission. Fig. 
H] indicates that the strongest boundary field is at r — r max and Fig. 3 gives 
f.o.m. ~ 1/8. 

Another relevant quantity is the group velocity t> g = cLu/dP, which can be 
easily calculated from (|14[) when the phase velocity w p h = lo/P is not imposed. 
For the u>a = 5.0 mode studied above, v g — 0.68. The guide can therefore be 
quickly filled with the wave and quickly emptied. To save energy, one may use 
a recirculated wave. 

We also have preliminary results for a waveguide with a stronger helical 
bending: e = 0.6, q = 1.5. Imposing v p h = 1 — 10~ 3 , three first modes are seen 
at ua = 2.2, 3.3 and 4.0 with \E z (0)/B z (0)\ = 0.35, 0.3 and 0.8 respectively. 
The number of sutture points, N ~ 10 — 20, was enough to obtain stable roots 
but not yet stable field amplitudes. This is like in the variational method in 
quantum mechanics, where the energy converges faster than the wave function. 



9 




-0.5 0.5 1.0 -0.5 0.5 1.0 



Figure 3: Cross section of the electromagnetic field J-(x, y) = F(0, x, y, 0) along 
the line y = for the TM-like mode at cua — 5.0 (circular base, e = 0.4 and 
q = 1). Left curves, from top to bottom : \£ r (x,y)\, \£ v (x, y)\ and \£ z (x,y)\ at 
y = 0. Right curves : same for B ri B<j> and B z . 
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Figure 4: 3- dimensional maps of |J5 r |, \E V \ (upper mapps), \E Z \ and \B Z \ (mower 
mapps) for the same mode as Fig. 3. 
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5 Conclusion 



This study has shown that a linear accelerator for ultrarelativistic particles may 
work with a helical wave guide. Until now we have made detailed numerical 
calculations only for a guide with circular base, one set of paramctres {e, q} and 
the lower TM-like mode. The figure of merit (accelerating field) / (strongest 
electric field at the boundary) of this mode was only about 1/8, but we hope 
that different guide parametres or different base shapes can give higher values. 

The study has also validated the truncation-and-sutture method. This eco- 
nomical method seems to be well suited to calculate the modes of a scalar or 
vector field in a 2- or 3-dimcnsional homogeneous cavity of arbitrary, but smooth 
enough shape. 



6 APPENDIX : The M matrix 



Defining 

C n j = coefficient of a; in the condition £ ■ a = at r„, 

Cn.i+N = coefficient of bi = ai + N in the condition £ ■ a = au point r n , 

C n +N,i = coefficient of a; in the condition £ ■ f = at r n , 

Cn+N.i+N = coefficient of 6; in the condition £ ■ f = at r„, 

the 2N x 2N matrix of the linear system of equations expressing the boundary 

conditions is decomposed in four N x N submatrices : 



/... 



M = 



C n .l 



V... 



+N,i 



l+N 



+ N,l+N 



...J 



C n .l 
Cnd+N 
C n +N,l 
C n +N,l + N 



[ipx-!a Sl + (p-u)) 2 ifj x +icr S; ] ( 



The matrix elements are, for I ^ 0, 

) 2 

[V'A-i^ + (p + ij x+l a -j e 
[si qr (tpx-i - (p - lu) 2 ipx+i) + 2(p - w) ip\ 
[si qr (tpx-i - (p + u) 2 tpx+i) + 2(p + uj) tp. 



illfr, 



and for I = 0, 



Cn,0 
C n ,N 
C n +N,0 
C n +N,N 



P ipi a r 
iuj ipi a v 

uoqripi . 
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In these formulas p = pi = P — Iq, k = ki = \uo 2 — (P — lq) 2 ] , A = |Z|, 
s; = sign(l), ipx = Jx(kir n )/ki and we have introduced the unitary complex 



(34) 



(35) 



(36) 
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number 

cr r (r») +ioy(r n ) _ ^_ llf>n <r x (r n ) +i(7 y (r n ) 



cr n = 7-r, — ^ = e 



CT(r„)| ' |a(r„) 
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